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ABSTRACT 


Total heat transfer coefficients were determined for a heated, 
copper cylinder placed horizontally ina fluidized bed. The fluidized 
particles were glass beads ranging in size from 215 um to 3.40 mm. 
Fluidization was accomplished using air superficial velocities from 
0.045 kg/me:S Gere. 555 kg/m@:S. Temperatures were measured using 
thermocouples and the resulting h values tabulated and plotted as a 
function of superficial velocity and relative velocity. Nusselt 
numbers based on particle diameter and tube wall diameter were also 
determined. The results show consistent patterns with other workers 


and are substantiated by empirical correlations. 
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INTRODUCTION 


Fluidized Bed Combustion is currently being studied intensely due 
to its very promising potential for producing clean and relatively 
inexpensive power from our nation's vast coal resources. Beds 
containing limestone or dolomite will absorb sulfur dioxide, thereby 
providing a sensible way to burn high-~sulfur coals to generate power 
while maintaining high air quality standards without the use of 
expensive wet scrubbing equipment. Furthermore, the long-term goal of 
a pressurized fluidized-bed coal combustor which will generate steam 
and compressed hot gases to drive both steam and gas turbines ina 
combined power cycle holds the great promise of increasing the plant 
thermal efficiency by fully several percentage points over 
conventional power cycles. 

Heat transfer between the combusting bed and the working fluid 
ordinarily takes place by means of submerged horizontal boiler tubes. 
Knowledge of the magnitudes of the heat transfer coefficients and 
temperature distributions along the tube surfaces is essential for an 
optimum design. The rate of heat transfer to a submerged tube depends 
upon a number of complex factors, including the properties of the bed 
material and the fluidizing gas, bed and tube geometries, and the 
fluidized state [1]. Many investigators have reported measurements of 
overall heat transfer between fluidized beds and horizontal tubes for 
a wide range of operating parameters [2-8]. 

Using a model of bed~-to-surface heat transfer in which the total 
heat transfer consists of the additive contributions of gas convection 


and particle circulation, the heat transfer behavior may be determined. 





Sance particle circulation iS difficult to measure, its contribution 
can be inferred from knowledge of the gas convective part and the 
total heat transfer coefficient. The gas convective component was 
measured by Matt Wade [9] in a separate set of experiments using a 
mass transfer analogy technique. The technique employs the use of a 
sublimating surface to indicate the convective component. The total 
heat transfer coefficients are determined in this paper and will be 
used for comparison to the gas convective results of [9]. 

Briefly, the investigation involves placement of a hollow, copper 
cylinder horizontally in a fluidized medium of glass beads. The 
cylinder represents a typical boiler tube and the glass beads are 
analogous to the boiler bed. Bed sizes were selected to model typical 
coal particle sizes. The cylinder has a heating cartridge placed in 
the hollowed portion which applies a uniform heat flux. Thermocouples 
are used to measure the bed temperature and resulting cylinder 
temperature. Knowing the heat applied and dimensions of the apparatus 
allow the total heat transfer coefficient to be calculated. 

The results of these heat transfer experiments will be compared 
to the work of several previous investigators [2,4,7]. In addition, 
the present results will be utilized in determining the ratio of the 


gas convective component of heat transfer to the total heat transfer. 





THE EXPERIMENTAL APPARATUS 


General Arrangement 


The experimental apparatus consists of a compressed air system, 
the containment bed, fluidized particles, the test cylinder and the 
instrumentation package. The air system produces, distributes, 
regulates, monitors and contains the working fluid. A general 
arrangement schematic is presented in Figure 1. Briefly, the figure 
depicts the oil-free compressor as the air source. This compressor 
delivers nearly 350 SCFM at 100 psia. This maximum flow rate was 
sufficient to provide fluidization over the entire range of particle 
sizes used in the experiments. 

The distribution system is composed of 2 in. PVC pipe and 
fittings arranged in a convenient air-tight geometry. Pressure 
regulators, pressure gauges and overall flow meters were used to 
monitor and control the flow. The air then flows through one of two 
rotameters for flow rate measuring. The rotameter selection depends 
on the particle size and associated fluidizing velocity under study. 

Flow control is achieved through the use of gate valves forward 
of the rotameter. The gate valve was selected instead of a globe 
valve to provide greater sensitivity in controlling the flow rate. 
After leaving the rotameter, the flow enters the diffuser where it is 
expanded to the dimensions of the bed itself. The pressure drop 
across the fluidized bed is observed by examining a mercury manometer 
tapped into the entrance of the diffuser. An air distributor system 


spreads the flow evenly to provide uniform fluidization. After 





passing through the beads, the air is exhausted to the outside by 


means of galvanized duct work. 


Fluidized Bed 


Fluidization is accomplished in an 11-inehn inside diameter 
plexiglass containment cylinder mounted vertically on supports. 
Figure 2 shows the bed in detail. The bed has a 6.5 inch x 9 inch 
removable access window with a matching leak-proof cover. The window 
allows the bed to be filled with particles as well as access to the 
test Section. 

The bead fill level is also indicated. This level is a constant 
regardless of the bead size under study. Two drain spouts allow the 
particles to be removed to independent levels. The upper level allows 
clearing only the test section while the lower spout is used to clear 
the entire containment of particles. The spouts are sealed with duct 
tape attached flush to the bed wall. This taping method minimizes 
flow disruption. The flow is left undisturbed by using smooth, 
polished inner walls and by having any protrusions away from the test 
section. 

Hnesbase or the bed has a conical, 2 inch to 11 ineh sheet metal 
diffuser acting to expand the compressed air. The length provides 
gradual, turbulent-free expansion. A distributor plate made up of a 
packed bed of 3.40 mm diameter glass beads sandwiched between two fine 
mesh brass screens forms the air distributor between the bed and the 
diffuser. The screens are 3 inches apart and secured by tack welds. 


The mesh of the screen is small enough to prevent the smallest 





particle size, P-010 from falling through it into the air distribution 
system below. 

The interface between the diffuser/upper screen and the bed is 
made by a gasket-filled flange connection. The gasket is a composite 
of 1/16 inch rubber material layers and tube type sealants. The 
flanged gasket seal maintains the air-tight integrity of the bed. The 
actual bed is shown in Photo 1. The advantage of using a plexiglas 
containment is apparent from the photograph. Namely, the extent of 
fluidization and general fluid behavior is clearly visible. Photo 2 


also shows the bed when fluidized. 


Fluid Particles and Rotameters 


The fluidized particles are industrial quality glas spheres of 
various sizes. The glass is made inert from its lead-free, soda-lime 
Silica composition. The spheres have uniform roundness and each size 
falls within a narrow band of dimensional variation. Since larger 
particles require successively greater volume flow rates to fluidize 
them, rotameters of different flow ranges must be used to maintain 
accuracy. Three meters were used having ranges of 1.35 SCFM to 27 
SCFM, 12.82 SCFM to 160 SCFM and 100.5 SCFM to 1005 SCFM. They will 
be identified as small, large and blue, respectively. The first two 
occupy interchangeable positions and are shown in Photo 3 and Photo 4. 
The blue meter is identified by its color and is the monitoring meter 
in the distribution system. Photo 5 shows the blue meter. This meter 


is used because its range exceeds that of the large meter and because 





meters of the same type as the large and small ones are not available 
in the appropriate size. 
Table 1 summarized the bead identifier, average particle size and 


the meter used. 


TABLE 1: Summary of Particle Name, Bead Sizes and Meter Used to 
Measure Flow 


BEAD 
EDENTIE TER 


AVERAGE 
SY 













2.15x1074 
5.125x1074 


9.45x1074 


2.03x1073 


3.40x1073 


The readings from the three meters were all converted to volume flow 
rates. In this way, the scale differences of the meters were 
eliminated and comparisons were made on a uniform basis. To better 
visualize the particle sizes, P-010 is about the size of table salt, 
P-O470 resembles large grain sand and A-340 is approximately the size 


of "BB's" for a pellet gun. 


The Test Cylinder 


A one~ineh diameter copper cylinder is the focus of the 


experiment. Its seven inch length and 3/16 inch wall thickness were 





designed to meet the heat transfer requirements of the experiment. 
Figure 3 and Figure 4 show the cylinder's design in detail. Copper 
was selected because of its high thermal conductivity, thus providing 
a uniform outer wall temperature when heated. To measure the 
temperature achieved and to verify the absence of circumferential 
differences, four chromel-alumel thermocouples were attached at 90° 
positions around the circumference. The thermocouple beads are 
attached below the wall surface and the wires run through machined 
grooves. The grooves are filled with thermal epoxy and sanded flush 
to the surface. In this manner, not only are the temperature 
characteristics unaltered but also the flow around the cylinder is 
left undisturbed. 

An electrical heater cartridge is placed inside the cylinder's 
hollow core. The cartridge acts as a radially symmetric heat source 
to generate the desired temperature gradient between the test cylinder 
and the bed. The cartridge matches the cylinder in both length and 
diameter. To ensure effective conductive contact between the 
cartridge and the cylinder, thermally conductive paste was applied to 
the cartridge before inserting it into the cylinder. The combined 
effects of the planned tight fit and the paste minimized any possible 
@Cenvact resistance. 

The cylinder is supported on two nylon end caps fastened into 
plexiglass clamps on the bed wall. Figure 5 shows the end caps or 
spacers. Paramount in their design is the matching with the copper 
cylinder of outer diameters, the fastening with recessed set screws, 
their length and the material selection. In this manner, the fluid 


flow is undisturbed and the test cylinder is positioned in the bed's 











center. Heat transfer end effects may be ignored because of the 
extremely low thermal conductivity of the nylon material. As a 
further precaution, the thermocouple wires and cartridge power card 
run through a hollow spacer and exit the bed directly. 

The assembled unit is shown in Photos 6 and 7. Photos 6 and 7 


show the cylinder removed and partially disassembled for clarity. 


Instrumentation Package 


The instrument package is relatively simple. It consists of a 
variac power supply, a radial, multi-terminal switchbox, a digital 
voltmeter, chromel-alumel thermocouples and an insulated icebath. The 
circuit diagram depicting the component relationships is shown in 
Figure 6. The actual equipment used is presented on Photo 8. 

The theory of operation combines the instruments to extract the 
temperature of the cylinder, the bed temperature and the power 
supplied. Other meaningful values may then be refined from the data. 
Briefly, the variac provides an adjustable current source to the 
heater cartridge. This source is indirectly measured by using the 
voltmeter. The voltage coupled with calibration information about the 
cartridge's resistance produces the power supplied to warm the 
cylinder. The digital voltmeter is a convenient and accurate 
measuring device for both AC and DC voltages present in the experiment. 
The DC potential differences arise from the use of thermocouples. The 
thermocouples use .010 m chromium wire and .010 m aluminum wire with 
teflon insulation. The thermocouples were manufactured in an inert 


gas environment by electric arc fusing. Manufacture using this 





technique assures measured values corresponding closely to 
standardized Type K temperature conversion tables. 

The radial switchbox is at the heart of the instrument package. 
It provides a common junction for all the thermocouples as well as 
providing easy selection of voltmeter inputs. Insulation packed in 
the switchbox minimizes false DC readings due to ambient temperature 
fluctuations. An insulated icebath is used as a reference for 
thermocouple voltage to temperature conversion. Careful use of 
distilled water and ice provided uniform reference from one data run 
to another. In addition to the cylinder and icebath thermocouples, 
three thermocouples indicate the inside bed temperatures. The 
thermocouples penetrate the bed wall at positions of 1/3 the bed 
circumference and 1/3 the elevation to the test section. This helical 
placement pattern averages out local temperature gradients to yielda 
good overall impression of bed temperature. 

Data collection using the instruments discussed requires manual 
recording of voltmeter readings. Automated "Dataloggers" on the other 
hand, electronically scan, average and convert the readings to 
temperature. Attempts were made using a Monitor Labs 9302 datalogger 
but the results proved unsuccessful. The reason for the failure is 
attributed to excessive static discharge. The static is generated 
from the plexiglass bed's natural interaction with the glass beads. 
The charge then travels down the thermocouple wires into the 
datalogger causing problems ranging in severity from erroneous 
readings to complete microprocessor failure. Mitigating efforts 
included: grounding the bed, grounding the datalogger, installing 


terminal connection slunts, electrically insulating the test cylinder 





10 


with 4 coats of polyurethane and even applying commercial static 
minimizer sprays. Failure of all attempts to reduce the static to 
within the datalogger's acceptable limits forced the abandonment of 
automated data collection. The digital voltmeter proved to be less 


sensitive to the static and was selected as the logical alternative. 
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FIGURE 1: Schematic of Overall Flow and Fluidizead 
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FIGURE 2: Detailed Schematic of the Fluidized Bed 12 
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FIGURE 6: Circuit Diagram for the Instrument Package 
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PHOTO 2: Close-up of Fluidized Bed Containment 
Vessel (Test Section) 
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THE EXPERIMENTS 


Procedure and Data Collection 


In an effort to minimize human errors an established experimental 
procedure was used. The standardized sequence steps were: 
1. Start compressor and wait until operating pressure of 100 psi 
was obtained (15 min.) 
la. Prepare icebath 
2. Adjust air flow rate to achieve desired fluidization 
3. Adjust variac to supply amperage where a temperature 
difference of about 10°C could be achieved 
4, Wait until a steady-state temperature is achieved (10 min.) 
4a. Check steps 2 and 3 
5. Collect data 
6. Repeat steps 2-5 
Data collection was also performed using a repeated series of events. 
To facilitate data collection and subsequent data reduction, a form 
depicting needed information was used. The form is presented in 
Figure 7. The form indicates locations to record voltage readings for 
the cylinder, bed, icebath and heater cartridge. The dial number is 
the corresponding terminal station on the radial switchbox. It is 
important to note the physical separation of the terminal stations. 
Separation prevents confusion and eliminates the possibility of 
erroneously recording data. The heater cartridge requiring AC volts 
to be read is separated from the thermocouples which require DC volts. 


The thermocouples are separated further into groups pertaining to what 





2 


temperature they are indicating. A reading is taken by selecting the 
dial number with appropriate voltage mode and observing the voltage 
displayed. This value is then recorded on the form. Fluctuations 
frequently occur because of the violent actions of the fluidized beads. 
When fluctuations occur, the changing values are observed and mentally 
averaged in an effort to isolate the true value. Errors of mental 


averaging are minimized by repeating the readings for five trials. 


Data Reduction 


For a given particle size at each flow rate, a separate data form 
Was prepared. Examining Figure 7 reveals locations to reduce the data 
into meaningful information governing the heated cylinder's behavior 
in the fluid. Reduction begins in the upper section by numerically 
averaging the readings at a given dial number for all five trials. 
The subsequent results are then averaged for the cylinder and bed, 
respectively. Using the icebath as a reference, the final bed and 
cylinder voltage values are calculated. This calculated value is then 
converted to a temperature by using the NBS Type K conversion table 
and linear interpolation. The temperature gradient is found by 
finding the difference between the cylinder and bed temperature. 

The heat supplied to produce the temperature gradient is found 


using the electrical expression. 


ye 
Rn 
Q = heat supplied 


Q = 


where: V = measured cartridge voltage 
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R = known resistance (53.45 ohms) 


Using the determined temperature gradient, heat supplied and knowing 
the surface area of the cylinder; the desired heat transfer 


coefficient may be calculated. The governing relation is: 


Qa 


fe 


where: h = overall heat transfer coefficient 


A = cylinder surface area (0.0142 m@) 
AT = temperature gradient (as measured) 
Q = heat supplied 


Nondimensional information is obtained by further data reduction to 
produce the Nusselt number based on particle diameter. 


The Nusselt number based on particle diameter is determined by: 


hd 


where: NUp # Nusselt number based on particle 


diameter 
h = overall heat transfer coefficient 
dy = particle diameter inom 
K = thermal conductivity of copper 


(2.63 x 1072 W/m: °C) 


The lower section of Figure 7 lends itself to programmed calculation. 


The large volume of data virtually necessitates programming. 
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To eliminate scale differences of the Nusselt number produced by 
using the various particle diameters as the characteristic dimension, 
the Nusselt number based on the outer diameter of the tube was 
determined. In this manner, a true indication of the total heat 
transfer coefficient is revealed for each particle size examined. The 
constant dimension allows direct comparison of Nusselt numbers between 
data sets. The Nusselt number based on the tube wall diameter was 
Calculated by multiplying the existing Nusselt number based on 
particle diameter by the ratio of the wall tube diameter to the 


particle diameter. The governing relation is: 


NUyt = NUD (=) 


D 
dp 


Nusselt number based on tube 


= 
> 
© 
e 
49) 
= 
& 

= 

ct 
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diameter 
NU, = Nusselt number based on particle 


diameter for each reading 


Dr = outer wall diameter of the tube 
(0.0254 m) 

dp = individual particle diameter in 
meters 


Air flow rates were standardized by determining the mass flow 
rate per minute through the fluidized area. This flow rate, called 
ae 1s in units kg/me:S, Flow readings from the three meters were 
converted to the appropriate units using procedures particular to the 


meter used to measure the flow. 
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All conversion processes involve knowing the crossectional area 
of the containment vessel. The area is determined by using the well 
known relationship for a circle. Based on a diameter of 11 inches and 
converting the result to me , the area is 0.063 m¢. This value becomes 
a constant for all calculations to determine "G". 

The flow rate for all readings taken on the Small Meter is 


determined with the following relationship: 


g - 2CeM 
A 
where: % = the decimal equivalent of flow read 

directly off the meter 

C = capacity of the meter (27 ft3/min.) 

0 = air density at atmospheric pressure 
and room temperature (1.201 kg/m3) 

A = crossectional area of the containment 
vessel (0.0254 m) 

M - metric conversion factor of ft3/min 


to m3/s (0.5886) 


The process governing standardizing readings from the Large Meter is a 
two-step method. The first step involves using a prepared conversion 
curve which permits the rapid conversion of the meter's units into 
SCFM. Briefly, the curve is a plot of percent of flow read off the 
meter versus SCFM. The conversion curve is presented in Figure 8. 

The figure is used by entering with the % of flow read from the meter 
during experimentation on the x axis and recording the corresponding 


flow velocity in SCFM off the y axis. This value is then used in the 





second step which is a direct 


performed using the relation: 


where: 


27 


calculation. The calculation is 


value determined from Figure 8 in 
SCFM 

air density at atmospheric pressure 
and room temperature (1.201 kg/m3) 
crossectional area of the containment 
vessel (0.0254 m) 

metric conversion factor of ft3/min 


to m3/S (0.5886) 


Flow rate determination for the Blue Meter is a simple calculation 


involving: 


where: 


b 


G 


2C pM 

A 

decimal equivalent of flow read 
directly off the meter 

capacity of the meter (1005 ft3/min) 
air density at atmospheric pressure 
and room temperature (1.201 kg/m3) 
crossectional area of the containment 
vessel (0.0254 m) 
metrie conversion factor of ft3/min 


to m3/S (0.5886) 
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Using the above procedures results in flow measurements having the 
uniform dimensions of kg/m@:S. 

The ratio G/Gyj;jn forms the dimensionless quantity which is used 
to indicate the relative magnitude of the flow rates used. The Gnin 
value is the smallest flow rate observed to fluidize the specific 
particle being studied. It then becomes the reference point from 
which to measure all other flow rates for that particle. Therefore, a 
G/Gnin ratio of 3 indicates the flow rate being used is three times 


the flow rate at the onset of fluidization. 





FIGURE 7: Standardized Data Collection and Reduction Form 
a eon 29 
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RESULTS AND DISCUSSION 


Results 


Data was collected for each particle size at numerous flow rates. 
The assembled results relating the: air flow rate (G lbm/min:ft°), 
the ratio of a chosen flow rate to the minimum fluidizing flow rate 
(G/Gain), the corresponding overall heat transfer coefficient (h 
W/m ° °C) and the dimensionless Nusselt number based on particle 
diameter (NUp) and Nusselt number based on wall tube diameter (NUyt) 
are presented by particle size. Tables 2-6 present the information 
for particle sizes P-010, P-0230, P-0470, A-205, A-340, respectively. 
Individual result summaries allow trends to be detected for each 
particle size. Relationships which link the results together are 
readily noticed from graphical representations. Figures 9, 10 and 11 
plot results for all particle sizes together. Figure 9 presents the 
relationship of h to G/Gnine Figure 10 presents Nup as a function of 
G/Gnine Figure 11 shows the relationship of h to G. Figures 13-16 
plot NUy,~ vs. G by individual particle sizes. Also indicated on each 
graph, are the predicted results obtained from published correlations. 
The Nusselt number based on tube wall diameter is graphically compiled 


on Figure 12. 
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TABLE 2: Results Summarized for Particle Size P-010 


Gnin = 0-045 kg/m@:S 

d = 2.15 x 1074 m 

k = 2.63 x 1073 W/mec 
G7G rah h NUp NUwe 
1.000 130.032 1.063 125.582 
1.388 193.741 14584 Se 133 
1.665 248 .656 25033 240.178 
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Results Summary for Particle Size 


Gmin 
d 


K 


G/Gmin 


1.000 
eees 
i570 


lee 


1.750 


e629 
1.848 
1992 
2.068 
2.202 
23h 
en (02 


3.068 


uoeuiwW 


0.256 kg/me-S 
5.125 x 1074 a 
2.63 x 107° W/mec 


79. 


168 


219. 


266 


Clic 


Zi 3- 


ete 


280 


aoe 


262. 


aie 


es. 


215 


680 


~ 347 


596 


se tO 


246 


B22 


oO 


O41 


Sf3 


021 


"29D 


368 


- 130 


NUp 


3 
12311 
279 
see 
OUD 
2320 
413 
SH 
444 
~ 496 
© 370 
poe 
ot 


p=0230 


NUwt 


TO 
1628 
alae 
eo 
202. 
263- 
268 
270. 
269 
alee 
266 
264. 


266 


968 
610 
O71 
122 
921 


962 


274 


454 


70:10 


387 


= b42 


011 


ee) 


33 





TABLE 4: 

Gmin = 

d = 

k = 

G G/Gmin 

0.564 1.000 
0.694 1.204 
0.774 1.344 
OmocT 1.435 
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Results Summary for Particle Size P-O47 
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Summary for Particle Size A-205 
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Results Summary for Particle Size A~340 
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Discussion 


Examination of the results supports the achievement of the 
experimental goal, that is, to identify the total heat transfer 
behavior for all particle sizes specific to the bed under study. 
Published values or correlations vary according to the physical 
make-up of the particular bed studied. Physical factors affecting the 
heat transfer properties include: size, type and weave of the 
distributor plate, the size, shape and thermal properties of the 
fluidized particles, moisture content of the air and lastly, the 
dimensions of the bed. The relationship between particle sizes and 
the trends of a given bead size are most readily noticed from 
graphical displays. Figure 11 shows the behavior of the total heat 
transfer coefficient as the absolute flow increases for all bead sizes. 
Figure 10 depicts the characteristic Nusselt number based on particle 
diameter for flow rates above minimum fluidization. Figure 9 shows h 
as a function of relative flow. All particle sizes are plotted 
together. To identify the individual particle sizes, different data 
symbols are used. The symbols are presented in Table 7 and on each 
graph for clarity. 

Examining the plot of h vs. G/Gpjn or Figure 9 produces 
interesting observations. The smallest particle sizes tend to have 
the highest total heat transfer coefficient. This behavior is most 
evident at G/Gmin values above 1.4. So, as particle size increases, 
total heat transfer decreases, in agreement with many previous 
investigations [2,4,7]. The decrease of total heat transfer with 


particle size is due to the increasing importance of gas convection 
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for the larger particle sizes. By tracing each particle size 
individually, a repeated pattern is noticed. The pattern tends to be 
one of rapid increase of h for a reasonably small change in G/Gpin 
until a maximum h is achieved at a specific G/Gnpj,. After this unique 
point, the h value decreases slowly for higher flow rates. Particles 
P-O47 and A-205 demonstrate the behavior most dramatically. Bead size 
P-010 was not observed in a flow rate range large enough to include 


its maximum. 


TABLE 7: Particle Name and Corresponding Data Symbol 






BEAD 
IDENTIFIER 








SYMBOL 


The uniqueness of a particular flow rate corresponding to a 
maximum h value is clearly seen in Figure 11. On this figure, the G 
at the largest total heat transfer coefficient may be read directly 
from the graph. The value of fluidizing velocity at which the maximum 
in h occurs is called the optimum velocity, Upp_. It has been 
observed by many investigators [10] and a correlation provided by 
Toades [11] may be used to predict its value. The correlation 
presents the optimum Reynolds number. The correlation yields the 


corresponding Uppt value to within a multiplied constant involving the 





a0 


specific kinematic fluid viscosity and the characteristic length of 


the apparatus used. The expression becomes: 


Ar rd 
Vopt = 78+ (5) (2avar) T) 


where: Ar = the Archimeder number 
Y = kinematic viscosity of air 
(1.55 x 109 m¢/S) 
L = characteristic length or the diameter 


of each particle size 
This relationship was applied for all particle sizes and the results 
compared to the experimental results of the present study. Table 8 


compiles the information by particle size. 


TABLE 8: Summary of Uppt by Particle Size 


Not Applicable 0.3580 


0.4225 One (2 
0.6886 0.8379 
licges00 leeeg5 


1.6919 1.5941 





The correlation velocities fit the experimental results very well for 


the largest two particle sizes, however the agreement is not quite as 





40 


good for the intermediate particle sizes. Total heat transfer 
decreases at velocities above the optimum velocity because of the 
increasing importance of gas convection heat transfer. As the 
fluidizing velocity is increased, a corresponding increase in bed 
porosity occurs. The excess fluid separates the particles which 
diminishes the strong contribution of particle heat conduction, 
thereby lessening the overall heat transfer. 

Closer examination of the h vs. G plot shown in Figure 11 
indicates that the dependency of h on G is stronger for smaller 
particle sizes. That is, the slope of the curve is greater for small 
beads. Using Figure 9, the steep rise is obvious when beads P-023 and 
P-O47 are compared to A-205 and A-340. For G/Gaj, values from 1.0 to 
about 1.6, the h values of the smaller particles changed from about 80 
to 265 and 135 to 236, respectively as compared to 150 to 195 and 165 
eo 190, respectively for the larger particles. 

This behavior may be explained by considering the mechanisms 
involved which compare the overall heat transfer coefficient. It is 
the sum of a conductive part and gas convective part. For smaller 
particles, the h value is dominated by particle conduction. Heat 
transfer predominantly caused the physical touching of packets of 
particles to the heated cylinder. Whereas for the largest particles, 
the overpowering mechanism is gas convection. 

Figure 10 plots the Nusselt number based on particle diameter 
versus relative flow rate for all particles. The plot clearly 
indicates that each bead size has a characteristic range of Nusselt 
number. The band tends to be narrow and at successively higher values 


as bead size increases. The increase in Nusselt number with 
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increasing particle size does not indicate conflicting results as 
presented in Figure 11. The Nusselt number calculation is controlled 
by the particle diameter dimension. This dimension overpowers the h 
value. The general shape of the curves for both figures is consistent 
with the earlier work of others [2,4,7]. 

Nusselt numbers based on particle diameters can be somewhat 
misleading. They may erroneously imply that heat transfer increases 
with increasing bead size. Nusselt numbers based on the tube wall 
diameter accurately portray the physical occurrence since the tube 
wall diameter is constant. Compare the implied behavior of Figure 12 
for NUyz as a function of G using all particle sizes to that of Figure 
10. As expected, the plotted curves have shapes identical to the 
plots of h on G. Figure 12 shows that the total heat transfer 
coefficient increases with decreasing bead size. 

Detailed examination of NUy-t in relation to G for individual 
particles is presented in Figures 13-17. Each display also has the 
behavior predicted by the three most applicable correlations. The 
correlations proposed by Grewal and Saxena [7], Andeen and Glicksman 
[4], and Vreendenberg [2] are shown below. Each investigator 
attempted to account for the many factors which affect heat transfer 


in a fluidized bed. 
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Grewal & Saxena: 


2 OND 


GD Co. Dnd/2g1/2,9-23 
ieee — stele UW PS ee 0.30 
wt 7( e)( Oru d 0 z) ( ke } Pr 
where 


= 4+ {4 
Pa aa [0 ‘Feels = TaeEID 


Andeen & Glicksman: 


GD ne 0.326 ne 
NU = 9OO0|1 - ° 
we = 900(1 ~ c)[ (PHB) (ee) br 


where: 
173 
1 200 G 
Som: {0.4 A ; } 
: oO los pr )g 
Vreendenberg: 
OFs 
Gap uc 
NU og (eS 0.3 
wie [ ( a Masyrg)! Rr 


The nomenclature common to all the correlations and the values 


which correspond to this experimental setup are as follows: 


E = bulk porosity (as determined above) 

G = superficial mass fluidizing velocity (G in. kg/m@:S) 

Dr = outside tube diameter (0.0254 m) 

Qs = density of the bed material ~- for glass beads (2500 Kg/m>) 


oF = density of the fluidizing air (1.201 kg/m3) 





43 


uo = fluidizing air viscosity (184 x 1077 N-s/m2@ 
dp = particle diameter as shown on TABLE 1 
g = acceleration due to gravity (9.81 m/s°) 
Cos = specific heat of the bed material - for glass beads (750 J/kg-°K) 
Ke = thermal conductivity of the fluidizing air (26.3x107° W/m: °K) 
eee = Prandtl number of air (0.707) 
ds = sphericity of the beads (assumed to be 1) 


For the smallest particle size, P-010, the Andeen and Glicksman 
correlation most accurately represents the data. Figure 13 shows the 
experimental NU, values and the correlation values as a function of G. 
This representation is about 10% different from measured values. 
Figure 14 presents the same information for P-023. Here, the Grewal 
and Saxena correlation coupled with the Andeen and Glicksman 
correlation form the upper and lower bounds, respectively. Both 
correlations could be interpreted to represent the data. The Grewal 
and Saxena correlation is most accurate for intermediate superficial 
velocities. At larger flow rates, the Andeen and Glicksman model is 
most applicable. It should be noted that none of the correlations 
account for the rapid increase of the heat transfer coefficient at the 
smallest flow rates for either P-010 and P-023. 

Figures 15, 16, 17 present the experimental and correlated 
behavior for bead sizes P-OH47, A~205, A-340, respectively. The 
experimental results are very poorly predicted by the three 
mathematical models in all cases. The correlations anticipate NUy¢_ 
values significantly lower than those actually achieved. The 


correlations were intended to be used for particles less than 500 
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micrometers. Since the last three bead sizes clearly exceed the 
particle size limit, extrapolation may not be valid. The results are 
qualitatively compatible and quantitatively substantiated by the work 


of others in the field. 
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FIGURE 10: : Nusselt Number Based on Particle Diameter 
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Tn FIGURE llt Heat Transfer Coefficient Against Absolute Mass Velocity 
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FIGURE 13: Comparison of Experimental to =mperical 
Nusselt Numbers Against Absolute Mass 
Velocity for P-010 
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FIGURE 14; Comparison of Experimental to Empirical 


Nusselt Numbers Against Absolute Mass 
VElOeciey for P-O25 
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PIGURE 15: Comparison of Experimental to 


Emperical Nusselt Numbers Against 
Absolute Mass Velocity for P-047 
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FIGURE 16: Comparison of Experimental to Emperical 
Nusselt Numbers Against Absolute 


Mass Velocity for A-205 Oo Ae os 


Grevual 4 
SOK EMNOAX 





__ Andee & 
Glickerra) 


rms 


- ~ 0+ eae ee a. | oe —__\lreancien~ 
wera 








es eh mya eeee pe —— 
oe eee eee — a - ee eee See ae ieee aes aay ae eae 
? 
~ -— ~ _ +r — i —~ _ 
j 
ated SS ee ey Ser STE eee ae a os — —- = ae 
t 
—_—— ~— — —_——  — eee oo _— — _ —_—— — — _ —— 
5 
| \ 
‘ ; 
i 
i 
‘ 
ao ee \© a rn ee ey Fe ps a een ge i a nae eee Se a es 
t 
» ' ’ 
\ , $ 
t i ; 
Soe ee ame aaa — = ee — — 
' 
ay | see | eennancons ee ae Waa ee ae Be x epee OT enn a ee — 
; 
ao? 
i eae — 
' onan 
= oe — = = — ae — a — 


—— ee — —- — i ee — = See tra ead as ks 
ee —_—- ee — ee eee —— ee ee te — ee 
‘ 
‘ 
_— — a — ——— — ee —— ee on —— = ———= ee 





, ’ 
a = eee oe et os —— — ae rot ore = 
So ee poe ne cena ~— — a —_ ee ee a 





= — 


FIGURE 17: 


_ 


ene ——— = a | See enn erro — 
ees Se A\@d er 
ae 


oe —— ed a aoe —~ 


—— oe - ee rene Sy 





ee ee ee LS ct ae ~_ 
~~ —— 


23 


ee ee ete -——— ee i oe — a 


Comparison of Experimental to Emperical 
Nusselt Numbers Against Absolute Mass 
Velocity for A-340 


- Eg cn 


ee ANNE RE snus Sane e peer een ee 
t ' ‘ . e 
I 


eat 


Grevoal & 
Sa kewa. 
-__ Avcleent 
Clicks man 


—\freertrcken 
Daca 








54 


CONCLUDING REMARKS 


Experiments were conducted to identify the heat transfer 
characteristics particular to the fluidized bed studied. The goal was 
to quantify the total heat transfer coefficient as a function of the 
mass velocity for five different fluidized particles. The particles 
ranged in size from 2.15 x 1074 m to 3.40 x 1073 m. Examination of 
tabulated and graphical data revealed the following generalizations. 

The smallest particle sizes tend to have the highest total heat 
transfer coefficient. So, as bead size increases, the overall heat 
transfer decreases. The decrease of total heat transfer with 
increasing particle size is attributed to the successively growing 
importance of gas convection. 

Each particle size exhibited a similar h versus G pattern. The 
pattern is characterized by a rapid increase inh for a small change 
in G until a maximum h is achieved at a specific G; followed by a 
Slowly diminishing h value over all remaining G values. The patterns 
first stage clearly demonstrates a stronger dependency of h on G for 
smaller fluidized bead sizes. The superficial velocity which 
corresponds to the maximum heat transfer coefficient in each case is 
the optimum velocity for that particle size. This optimum velocity, 
Uopt Was experimentally identified for each bead size. Knowing this 
ideal point will aid subsequent heat transfer experimentation using 
the same apparatus. The experimentally determined Vopt values compare 
favorably to empirically defined values found from the Todes 
correlation [11]. The last section of the pattern shows h slowly 


tailing off. The gradual decrease is caused by the increasing bed 





DD 


porosity at higher velocities. The larger fluid quantities decrease 
particle heat conduction, thereby enabling gas convection to dominate 
the heat transfer process. 

The dimensionless Nusselt number was used as a basis for analysis. 
The Nusselt number based on particle diameter, NUp#, shows that each 
bead size has a narrow range of NUp# values with the values increasing 
with particle size. The Nusselt number based on tube wall diameter, 
NUy-#, was used to directly compare the heat transfer data at all 
velocities among the various particle sizes. NU,,# behaved in a 
manner identical to the behavior of h. The validity of the NU,;+# 
values was tested against results obtained from numerical correlations 
by Grewal and Saxena, Andeen and Glicksman, and Vreendenberg. The 
correlations substantiated findings for the smallest three particles. 
For the largest two particles, the correlations do not reflect the 
effects of the large beads. This makes them invalid and were not 


applied in the analysis. 





nO 


Pa 


56 


REFERENCES 


Botterill, J.S.M., Fluid-Bed Heat Transfer, Academic Press, New 
Work. 1975 < 


Vreedenberg, H.A., "Heat Transfer Between a Fluidized Bed and a 
Horizontal Tube," Chemical Engineering Science, Vol. 9, 1958, 
Epa 52-00. 


Chen, J.C., "Heat Transfer to Tubes in Fluidized Beds," ASME 
Paper 76-HT-75, ASME-AIChE Heat Transfer Conference, St. Louis, 
MO., 1976. 


Andeen, B.R. and Glicksman, L.R., "Heat Transfer to Horizontal 
Tubes in Shallow Fluidized Beds," ASME Paper No. 76-HT-67, ASME- 


AIChE Heat Transfer Conference, St. Louis, MO., 1976. 


Krause, W.B. and Peters, A.R., "Heat Transfer Mechanisms Near 
Horizontal Heat Exchange Tubes in an Air Fluidized Bed of 
Uniformly Sized Glass Particles," ASME Paper No. 79-HT-88, ASME- 


AIChE National Heat Transfer Conference, San Diego, CA., 1979. 


Grewal, N.S. and Saxena, S.C., "Effect of Surface Roughness on 
Heat Transfer from Horizontal Immersed Tubes in a Fluidized Bed," 


ASME Journal of Heat Transfer, Vol. 101, 1979, pp. 297-405. 


Grewal, N.S. and Saxena, S.C., "Heat Transfer Between a Hori- 
zontal Tubes and a Gas-Solid Fluidized Bed," International 


Journal of Heat and Mass Transfer, Vol. 23, 1980, pp. 1505-1519. 


George, o.b. and Grace, J.R., “Heat Transfer to Horizontal Tubes 
in a Pilot-Seale Fluidized Bed," Preprint, ASME-AIChE 20th 


National Heat Transfer Conference, Milwaukee, WI., 1981. 


Wade, M.L., M.S. Paper, "Contribution of Gas Convection to Heat 
Transfer between a Submerged Horizontal Tube and a Fluidized 
Bed," Pennsylvania State University, 1984. 


Goroshko, V.D., Rozenbaum, R.B., and Todes, O.M., Izv. Vysshikh 
Ucheban. Za vedenii, Nef 4; Gas (Baku) No. 1, 125, 1958. 


Toades, O.M., In "Applications of Fluidized Beds in the Chemical 
Pidiucury, “Part If, pp. 4-27, Izd. "Znanie,“ Leningrad. 





- 











| Thesis 20 28 f 
tT) 64 terepe lat t i 
Overall heat transfer 
coefficients for a hori 
zontal cylinder in a 
fluidized bed. 


T2764 Terpolilli 
Overall heat translrer 


Poet lelents for a hori- 
Pet ieee | cee Dee 
fluidized bed. 





| lel Be 








be 
CANS ae rae i 4 ee 
yet Bee MON : 
es a pula hea thesT2764 
VE aie 


iH 
my Pi tf 
ey aia eae SSE FiGM eS f 4 a3 ind ee ee 
oe bate May és Bue a | {| HT 
Cae fe RELA RORY 1 pay 

- \ Um +4 i y Y' iP 

rp k aed i nan S OX f 

Lit | | 


eon 4 iia: a ey ig 


eh | ~e@eal 
Ce a de oe 
- ee hid oe Cra + 


> 


£4 aE Ty Deh 1 
pin: Wi en Seeme § Overall heat transfer coefficients for 


1 


a frye? Be ee re | Op See ee 4 Og 

atk BR a 4 988 A a Ld A 

ee ce | ae p Beanie 3 2768 002 03446 4 
yn get at 


ee 

oy it ae 

oy ri} road 
RO Ay Tas 

rises cian 


TNE g 
We 


es 
te a “obi 
i 


he raes oe 


an sk 

ae Cd 

aie ie tae Tw, ett ices Afar eet 

iat oe ne ae Ny a We ‘a rt ie fate rails 
ft} NE: iw is Pb wed 


a 


pony ' ire bad 

i en Pan pen SERRE eRe 
Coat Lease es te , 4 TY tne 
vial dep Piatt rrk ih f ie ay AP hanet oe 5 See 


i 4 
bana rs ere pe ’ Hoaa Sd Wie ye aed +i ea oar ar i. 
BI é i Skate KE i iti as iD Bre id rate TT Ria. re es 4 o “ ase nant H yon 


LE nN Cot io 
“ aya) 
aren DOM 





Fi oe om ” " 
Sy et Bik: ; ee se ae 
alee ya Ra ay pe? 
in tatie ar eke 

ava rH 


ets oe ey 

ysieones ee neath * 
SOT ena Lee Ae Lea : ier a WN te 

f te vn ¥ L * i yo a 0 

Ay the Boe ae i a a he ee 


ars 


Aor an ‘s a 

Eire cramer ony ol ert) 

oe oe Sir eae mek RRL rie 
a a Ta erry Pett ; wr pi 
a ee ine Ales 


i 
en 
‘a Poet 

ees 
ue tia me yt Pena erica at ea ar 
ra . ; 
Rate Mae Bn ae actit ate i Re 
ede oe %, LY: Mees Pert ae is ve 
x yay ys CaP Te 
ar Es Fe hintn hi AN iY ae nt a alas role 


by an on Me (ty 
iran Ra Ne meee 


e, Uv 
ve var tant 
ae beat eat Ps Hie 
rt i iM 
Y re ar 5 
ane ba - ey fa) 


oe ree, 
aon He rit 

ies i ee? 4 ne Poet te aoe i if a ty 
LOR han ee TL hae ort ret ke Se 
HY Rae RraL ls hy + ki ra is tv nae po ate at eae F 
es F “ he PLT oar * 

wy Bis tn if cs eg ae 

4 wie er ord ; 
ite Tat v ert i fe pas 


HIN Haiti 
ne hw A] Ui one 
Suge eet Hy , enn Py ee 















ai 
ww CL F) 
eo 


A iys™ WG Ht) Tae 4 Pi 
RUN rk pS Ria 































































































































































Setar 
Lape rh rae, 5 ‘ my tu 
Leu re a es Mead ea aa Broa ny 4 ye te 
De ae mM tia Ns tis ant in a Cai ae i ne oe whe , Hl 
pe) eats i ey ye aH 7) ye re bee o ro : 5 rr 
“H eats Uae ante cong oh YU ’ ak per Lor aes a che ed 
Mare pee ‘ Ls “a 
Kia by ks ead ra ae : & és ie Mioeac r u eel 
a a im er Sh tile it i 4 
Cet Sere i at ATH yet 5 eat) A aA so ‘ Pere oi 
sa SUR A Ae a ne iy ‘tte te oy yl i Mi ry Pra ry . We rs 5 
vi Seas ie it ‘i 4 ¥ ba Cue fe 
aN My Shae fea Dats , ee Af bas Te HF y AAT ae 
Sie) ts Cd NT acer ON ‘nt! ia Reet j ee e Pay , 
i A i Aes eee, ‘4 5 we a * 
| Ps eis, Or me: a id eae ah " st fii a 7 4 ° Goat he 48 
as ah ng rw a ce han ry oD p P oa ' ener ac | SE rn 
ua a i nes aH ae KY Sa yes he a ath A 
bg 4 : vt n 1p ee ihe - a fe 
: eens oR Aha , ete EAC APLC snr a Re taea ae 
a Ta Rae atk th So ES Wl aver oy : ae 
oa yet re toad Rea aed Ree aie Sir aa f on NE Ai rye qe Se wet y P 
ca. oe SN ae Han ig PAY ere cite a Pa bs i Pre . ee : n 
etre A PSN a ence ' 
NGA Hons, WaT) ae Sei? wa YY awe an ts 4 mm ; : 
ce MO kot) ss add i coy ei 4 a Set a MY f e +4 ne PI us 
rt Ld AL d Sa) 
ey Me ee SENS mat ea: 7 be Catia me way | a ae rt ~ihad ses 
ry rh ; , 7 
Pecos ty a ei aoe aha RAN e: nh m puetee ae ewes al 
Pod re é RN é, Td irs as 
~ baa Pee on a) g PF a Pu ‘ Ki 
‘eH oe Res Ate TS Pkt es ° + Da 
a4 Dee bi) rie Ye ee 6 vi ae Baie y Wap 
4 : Si is ee Le are aS 4 | _ ; 
h ead te hs p AL eee, baht ae te ial Fy By 
ry fs tet Sey Pity cote ty , ae 
BS ‘7 maT? ar el | 
oh te ad re he ot ? , “wy # 
is i! J La , a 
i , Pe a a 
i‘ Ce ah a se : oe ae 
ws Cet oe a re Pal at - "ih wat aee 
oN ean } ne sn Ng oe , "— _ F =! +b ae 
aaa . * 
mins ty oe ty tS Ny a ee Td , * 
i Ky Ly to Who yet am aa ie me 5 
‘ eo, “ P 
bas , * 2 ry ° 
Ed Ae Y | 
Nak wi Zig) a , Chk Aa! 5 Ft \ oe 7 5 Pa 
r hs SUSE HB ea B ac a a rx. % 
? a4 4 ) Le sy * a * rf L 
Tt % ra : ’ U 
a es ee M 
Xs ah Peg aa) Mi ae eg . 1 
~ Soe aay V8 pc ry a Th 4 i : 
ae 
vs ey ve se Sea Tid aan i ' 
agp, ot Py, eS , a 2) A al 
, iy beer Pan itad oth wih ny ig 
mn pina § EO a po one, 
aa no oe p erie ay Kho: o y PT A 
MN) ae Ps Ie Noe ae aaa ot al 
eS Aad a mA an ; Py " » : => cry 
; ory a) r 
ae | eee oe | : 
Re Hk ee ee ee 
wage. oF Aah eh a area 
a sy ie SR eg Pi Peer SS 
eG td rer 
A % ue eh ch’ a - OS Tey . 
te st a q pee aT a YP on us” ry a Me ; a ne LA ‘ 
hee , eh as AN So oF 
OCG taht A ut ra ar OR ‘ 
A recn bey on a 4 . My ‘ . 
haa ~ yey ie PUR NE oe | 
Se mi SF, . OEY + a ro 8 a a 
ave Pit > 4S Wh nureat oe ios un a 
ce ar Pi 
+n PRs aa af rf in t 
Pe ok hii Mh Y . 
A ees 
a » & «8 
Rta a 
Phy EEN i ry ‘ ore) in 
H athe ate + Real MOY ee \ Peli ) 3 
; Sans : wi BRE SAN Ce Pee ryt me Ane = 
SSeS rs wi MAD Oewe (h i > | Ni? Boh ios Rn hy be a N 
AN ” x Se . i aa a ca C ae ar ae i b 
G S D ‘i 7 ‘3 x Ra ans rt re OT ae] ats 's %, rs +E ' Peer 
. a > es I oa ut & ne Ye Ai ut a bs 
ASS enn ee i eae i . hh co eae uae baat : 4 ‘s ae x 
; ae as » 
Aare ; 8 Si a n y Saas A ‘ia . 
MZ a a Eka f eee Ohh Ca ha ” 
Se N uf ; Hay u ‘ 7 ys i wk ARN cm oats Nem f a) 
< \ : N > A ae Cl Nae ry, won Ss TCTs eG oe PO An Se % © fy hi 
‘ Hi oe bey Fiaon a or At ae ) is Pei a 
at SNe te iy ar i . 
a a FEO th Se Na re is é . - \ 9 by 
+ A ie, re vt Mane ia Rak ed os te rm a ‘al 
eh Tet aD a Sa ae p . 
BY Tayegnny Atta he 
By A eR te as .e a a 
Se 
ie Se Ae AS # ) A 2 b 
Fe pr nena ane , 
MaTe age, Face ee a. F 
ao ag rs ay ry bs “ . fs rt : wa 8 * 7 % as 
' \ b bs 
RATS, Sek Cat en a A. a 


DUDLEY KNOX LIBRARY 


re ore eb 4, 
wa, 
ty) 












Per 2 


